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ABSTRACT 

We discuss the effects of particle acceleration at the bow shocks expected in the binary pulsar system J0737- 
3039, due to the wind from pulsar A interacting with both the interstellar medium (ISM), and with the mag- 
netosphere of pulsar B. In this model, we find that the likeliest source for the X-rays observed by Chandra is 
the emission from the shocked wind of pulsar A as it interacts with the ISM. In this case, for favorable model 
parameter values, better statistics might help Chandra marginally resolve the source. A consequence of the 
model is a power law high energy spectrum extending up to < 60 keV, at a level of ~ 2 x 10~ 13 erg cm" 2 s" 1 . 
Subject headings: pulsars: general - pulsars: individual (PSR J0737-3039A,B) - radiation mechanisms: non- 
thermal - gamma-rays: observations - X-rays: binaries 
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1. INTRODUCTION 



The double radio puls ar system J0737-3039 tLvne et alJ 
2004; Kas pi et al~l 12004) is of great interest as a remark- 
able laboratory for probing strong field gravity and magne- 
tospheric interactions. It has also been detect ed in a 10 ks 
Chandra observation (McLaughlin et al. 2004), with an X- 
ray luminosity of Lx ~ 2 x 10 30 (ii/0.5 kpc) 2 erg s" 1 in the 
0.2- 10 keV range (where d is the distance to the source), 
and a reported X-ray photon number index of T = 2.9 ±0.4. 
The spin-down luminosity of pulsar A, which is expected to 
be channeled ma inly into its relativistic wind, is Ea — La — 
6 x 10 33 erg s" 1 (iLvne et alJEOOfl iKaspi et al. l l2004). Since 
L^~3x 10 3 Lx, only a small fraction of La is required in or- 
der to produce the observed X-ray emission. Since only 77 ± 9 
X-ray photons were detected, the determination of the spec- 
tral slope is difficult, and might be consistent with a flat vF v 
(r ~ 2), as expected from shock acceleration. Here, we ex- 
plore whether particle acceleration in the bow shocks of the 
pulsar A relativistic wind can explain the properties of the X- 
ray emission. The bow shock on the the magnetosphere of 
pulsar B involves only a small fraction of the pulsar A wind, 
due to the small solid angle it extends as seen from pulsar 
A. Therefore, it must have a very high radiative efficiency 
in order to explain the observed X-ray luminosity. On the 
other hand, the bow shock on the interstellar medium (ISM) 
involves most of the pulsar A wind and thus allows for a sig- 
nificantly smaller and more realistic radiative efficiency. We 
evaluate the expected high energy emission from this shock 
model, which also predicts emission up to tens of ke V. 

2. EMISSION FROM THE BOW SHOCK ON THE ISM 

At a sufficiently large distance from the double pulsar sys- 
tem, a bow shock forms due to the interaction of the wind 
from pulsar A with the interstellar medium (ISM). 2 This sit- 
uation is similar to that for a millisecond pulsar with a close 
low mass binary companion ( Arons & TavanD ll993l) . as far a 
the interaction between the pulsar wind and the ISM is con- 
cerned. The relative velocity of the center of mass of the bi- 
nary pulsar with respect to the ISM is 140.9 ± 6.2 km s" 1 on 

1 also Pennsylvania State University, Dpt. of Astronomy & Astrophysics, 
Dpt. of Physics and Center for Gravitational Physics 

2 The spin-down power of pulsar B is ~ 3 X 10 3 times smaller than that 
of pulsar A, so that its wind should have a negligible effect on the bow shock 
with the ISM. 



the plane of the sky (Rans om et all2004l) . A velocity compo- 
nent along our line of sight could lead to a larger total velocity, 
v ex t = 200v'2oo km s" 1 with v'200 ^ 1 ■ The head of the bow shock 
is at a distance R from pulsar A where the kinetic pressure of 

the wind balances the ram pressure of the ambient medium, 
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where p exl = n exl m p and n ext = hq cm are the ambient mass 
density and number density, respectively. 

Pulsar winds are thought to have a pair plasma composi- 
tion, perhaps with ions in restricted latitude sectors, and a high 
asymptotic bulk Lorentz factor (perhaps as high as ~ 1 4 — 1 6 
in the Crab Nebula and other young pulsar wind nebulae). For 
simplicity we assume a pure e pair plasma which holds a 
fraction e e w 1 of the internal energy behind the shock. We use 
a fiducial value of j w = 10 5 7 H ,.5 for the wind Lorentz factor just 
before the shock, however our main results are rather insensi- 
tive to the exact value of j w . We assume j w 3> 1 throughout 
this work. 

The ratio a of Poynting flux to kinetic energy in the wind 
is believed to be a ^> 1 at very small radii, while low val- 
ues of a <C 1 at large radii are inf erred from observations 
(e.g. (7-v3x 10" 3 for the Crab, iGaflant & Aronsl [l99i 
Spitkovskv & Arons 2004). It is hard to estimate the value of 
a at intermediate radii, which are relevant for our purposes. 
For the bow shock with the ISM which is at a relatively large 
radius, one might expect cr < 1. The shock jump conditions 
imply that the fraction e# of the internal energy behind the 
shock in the magnetic field is eg ~ a. However, amplification 
of the magnetic field in the shock itself could produce e# ~ 1 
even if a <C 1. Conversely, for a > 1 magnetic dissipation 
behind the shock might decrease the value of eg and make it 
close to unity. Therefore, we assume eg ~ 1, and to zeroth 
order we neglect the effect of the magnetic field on the shock 
jump conditions. 

In order to estimate the emission from the shocked wind, 
we will use the values of the hydrodynamical quantities at 
the head of the bow shock. To first order we neglect the 
orbital motion of pulsar A. The proper number density in 
the wind, as a function of the distance r from pulsar A, 
is n w = LA/47rr 2 m e c 3 7 2 .. The shock jump conditions at the 
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head of the bow shock imply that the shocked pulsar A 
wind just behind the shock moves away from the shock at 
j3 = 1 /3, and has a proper energy density e mt = La/2ttR 2 c ss 
1.3 x 10~ 9 noV2 00 erg cm" 3 and a proper number density n = 

2 3/2 7 w n„. = e int /(m e c 2 j w /y/2) = 2.3 x 10" 8 « V2 o7w,5 cm ~ 3 - 

This implies a magnetic field of B = 1.8 x iO~ 4 n 1 J 2 V2oo^J 2 G 
(in the fluid rest frame). The e ± pairs are assumed to be 
accelerated by the shock into a power law energy distri- 
bution dn/d r y e oc r y~ p , with 7„, < 7e < 7 max . Observations 
of synchrotron emission from electrons accelerated in rela- 
tivistic collisionless shocks typically imply p ~ 2-3. The 
average random Lorentz factor of the shocked electrons is 
(pf e ) = j w / \/2, and the minimal Lorentz factor is given by 3 
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where g = 3(p-2)/(p- 1) equals 1 forp = 2.5. The maximal 
Lorentz factor, from the requirement that the Larmor radius 
Ri = jeiiiuC 2 jeB does not exceed the width -qR of the layer of 
shocked fluid, is 



7max,l — 



eBr/R 



= 1.2 x 10 7 4 /2 (t?/0.3) 



(3) 



Here, the value of 77 can be estimated by equating the par- 
ticle injection rate into the hemisphere containing the head 
of the bow shock (9 < 90°), N/2 = L/2~f n m e c 2 , to the flow of 
shocked particles behind the shock outside of this hemisphere, 
2nr\R 2 nu, where n = 2 3 < /2 7lv n vv . = N/ V2nR 2 c and u = 7/3 are 
the proper density and 4-velocity (in the direction perpendic- 
ular to the shock) of the shocked wind at 8 = 90°. This gives 
r) w 1 /(2 3 / 2 m), so that 77 < 1 implies /3 > 1/3. At = 90° we 
expect > c s /c w 3" 1/2 and u > 2" 1 / 2 so that 77 < 1/2. On 
the other hand 77 < 0. 1 requires u > 5/ \pl « 3.5 which begins 
to be highly super sonic, and is therefore not very reasonable. 
Hence we expect 0.1 < 77 < 0.5 and use a fiducial value of 
77 = 0.3. 

The dominant emission mechanism is synchrotron radia- 
tion, and inverse Compton scattering can be neglected. The 
Lorentz factor of an electron which cools on the dynamical 

200 s, is given by 
4.7 x 10 10 

lc=-^T— = ^n • (4) 



time, f dyn ~ R/(c/3) = 4.9 x 10 5 n 
6irm e c 



<7t 



Bh 



1/2 

dyn e B n ' V200 



The synchrotron spectral break frequencies corresponding to 

7m , 7 C and 7 max are 
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We have vF v oc u 4 ' 3 for v < v m , vF v oc u®~ p " 2 for v m <v< 
mhn(v c , v mix ), and if i/ m - dx > v c (which is relevant for the next 
section) we have vF v oc v®~ p " 2 for v c <v < v m ^. 

Since r y c > 7 max , all electrons radiate only a small fraction of 
their energy. The fraction of energy radiated by an electron is 

3 More generally, this expression should be multiplied by a factor of (1 + 
Pp/p e ), which can be as high as m p /m e in the limit of a proton-electron 

plasma. Also, the factor of is valid for p > 2, while for p = 2 it should be 

replaced by l/ln(7max/7m) so that g = 3/ln(7 raa x/7m). 



~ min(l,?d yn /? c ) = min(l, 7e / 7c ), where t c = 6irm e c/ ' ajB 2 ^j e is 
the synchrotron cooling time. Averaging over the power law 
electron energy distribution, we obtain the total fraction e lat j 
of energy in electrons that is radiated away. For -f c < 7,,, (fast 
cooling), e ra d ps 1, since all electrons cool significantly within 
f d yn- For 7,,, < 7 e < 7 max , 
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For our fiducial parameters and p ps 2 we have e ra d ps 4.3 x 



10" 



1/2 

'0 v 200£ 



3 /2 

b (77/O.3). Most of the radiated energy will 
be emitted near f max at tens of keV. The fraction fx of the 
radiated energy in the 0.2- 10 keV Chandra range is ap- 
proximately given by the ratio of the average vF v value 
in the Chandra range (equal to the vF v value at some fre- 
quency vx within that range) to the peak vF v value. In our 



case, fx ~ (^/iw) (3 -" ,/2 « O.llnJ'v^ % 7 W0-3)- 1 
(this expression holds for vx < f m ax < v c ). There- 
fore, the ratio of the expected X-ray luminosity Lx = 

fxte^m&LA ~ 7 x lO^n'/V-joo 6 ^ 4 er S s_1 m tne Chandra 
range, to the observed L£ bs « 2 x 10 30 (fi? /0.5 kpc) 2 erg s" 1 is 

~ 0.35ny 4 V2ooe^' /4 (t//0.5 kpc)" 2 . This ratio would be unity, 
e.g., forno ~ 60 with v'200 ~ l,orfor«o~ lOand V2oo^2.5. If 
d « 1 kpc instead of sa 0.5 kpc, then we would need «o ~ 10 3 
and V200 ~ 4, which are less likely. Conversely, the con- 
straint is easier to satisfy if d < 0.5 kpc. According to this 
interpretation the high energy emission should peak at tens 
of keV (near is max that is given in Eq. with a flux of 

~ L x bs /f x 4nd 2 - 2.5 x 10- 13 n„ /4 v 2 ^4 /4 (y°- 3 ) erg cm" 2 s" 1 . 

Another contribution to the X-ray luminosity might be ex- 
pected from the shocked ISM in the bow shock. The energy 
injection rate is 4 ~ (v ext /c)LA ~ 10~ 3 La, which is of the or- 
der of the observed X-ray luminosity, and perhaps larger by 
a factor of a few. This could account for the observed X-ray 
luminosity, provided that fx^md ^ 0.1 -0.3. Here, the dy- 



namical time is fd yn ~ R/v ext and hv c = 22e B 3//2 n 1//2 V2oo keV, 
while v m is very low (in fact j m ~ 1), and the expression for 
t'max is the same as in Eq. Qwith the difference that here i^R 
is the width of the shocked ISM layer (instead of the shocked 
pulsar wind), and that es might be different (probably some- 
what smaller) in the shocked ISM. One might expect insta- 
bilities near the contact discontinuity between the shocked 
wind and the shocked ISM, both of the Rayleigh-Taylor and 
Kelvin-Helmholtz types, which could bring the magnetic field 
in the shocked ISM close to equipartition. For V200 ~ 1 and 
no > 10, vF v peaks in the Chandra range, so that we can have 

4 Comparing the energy injection rate per unit area into the wind ter- 
mination shock, L^/AitR 1 , and into the bow shock going into the ISM, 
(l/2)/9cxtv^ xt , and balancing the two ram pressures, La/4ttR 2 c and PcxtVj Xt , 
respectively. 
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fx ~ 1. Since the shock going into the ISM is Newtonian, 
one expects p«2, as in supernova remnants (SNRs). For 
«o ~ 60 this would imply e ra d ~ 0.2. From modeling of colli- 
sionless shocks in SNRs, which propagate into a similar ISM 
with similar shock velocities, a typical value of e e ~ 0. 1 might 
be adopted. The resulting value of Lx ~ (v ext / c)fxe e £radLA ~ 
2 x 10" 5 L A ~ 10 29 v 30 o(e f ./0.1)(/xe rad /0.2)( fl f/0.5kpc) 2 erg s" 1 
is only ~ 0.05L^ bs . Thus, this emission component may not 
easily account by itself for the Chandra observation (unless 
e e ~ 1), although it can contribute to that from the shocked 
wind of pulsar A. 

We note that Eq. ^ implies that the angular distance be- 
tween the double pulsar system and the head of the ISM bow 

shock is 6h s = 0.65(^/0.5 kpc)" 1 /!^ 1 ' 2 v^qq arcsec, and the rela- 
tively bright emission from the bow shock could extend over 
an angular scale a few times larger than this value. This an- 
gular scale may be resolved with Chandra, with longer in- 
tegration times, even though in the 10 ks Chandra detectio n 
it was reported as a point source (McLaughlin et al. 2004). 
If resolved, one might constrain the source angular size to 
< 1 arcsec. However, we note that the observed X-ray emis- 
sion is best explained from the bow shock with the ISM 

if n I//2 V2o ~ 0.35e^ 4 (of/0.5 kpc)" 2 , which in turn implies 

6bs ~ 0.23e 3 J 4 (d /0.5 kpc) -3 arcsec, that may be difficult to re- 
solve with Chandra unless d < 0.5 kpc. On the other hand, 
this suggests that d > 0.3 kpc, as otherwise the source should 
have already been resolved by Chandra, despite the poor pho- 
ton statistics in the current observation. 

The emission from the bow shock with the ISM is not ex- 
pected to show significant modulation at the spin period of 
pulsar A, Pa = 22.7 ms, or at the orbital period, P or b = 2.45 hr. 
The former is because Pa is ^6-7 orders of magnitude 
smaller than R/c. 5 The latter is because the orbital velocity of 
pulsar A is v OT b ~ 300 km s" 1 <C c, and the distance between 
pulsars A and B is R AB = 8.8 x 10 10 cm < R. 6 

3. EMISSION FROM THE BOW SHOCK AROUND PULSAR B 

Balancing the ram-pressure of the pulsar A wind with the 
magnetic pressure of pulsar B, assuming a predominantly 
dipole magn etic field, and a surface field strength of B» = 
1 .2 x 10 12 G (iLvne et all2004l) the distance of the head of the 
bow shock measured from pulsar B is R/, s «6x 10 9 cm. This 
is sa 0.07 of the separatio n betw een the two pulsars, Rab = 
8.8 x 10 10 cm (iLvneetalJ 12004 . Therefore, as seen from 
pulsar A, the fraction of the total solid angle subtended by 
the bow shock is VL/A-k = C^{R bs / Rab) 1 I ^ ~ 10" 3 C, where 
C ~ a few, its value depending on the exact shape of the bow 

5 This has two effects. First, any variability in the wind with the period P& 
will be strongly smoothed out by the time it reaches the bow shock. Second, 
the distance of the bow shock from pulsar A varies, with AS ~ R, so that the 
phase of the pulsar A wind that impinges upon it at any given time changes 
by ~ 10 6 — 10 7 periods. Since the same holds for the observed emission, it 
significantly averages out a possible modulation with a period of Pa , even if 
it exists in the local emission from a given location along the bow shock. 

6 The orbital motion of pulsar A affects the bow shock with the ISM 
mainly in two ways. First, the distance between pulsar A and the head 
of the bow shock changes by ~ ±Rab/2, changing the ram pressure by 
Ap/ p ~ 2Rab/R ~ 10~ 5 — KT 4 . Second, the wind is highly relativistic and 
behaves as radiation, so that its intensity in the bow shock rest frame scales 
as the fourth power of the Doppler factor 8 1 + /3 or |, cos 9, and will change 
in the range (1 ± /3 or b) 4 , resulting in a relative amplitude of ss 8/3 or b ~ 0.8%. 

Since (R/c)/P 0I t, 18n ' ^ v^oo, some additional averaging can occur due 
to the different phase of this modulation over the different parts of the bow 
shock, although this effect is not very large for our most promising model for 
which R is smaller by a factor of ~ 8 compared to its fiducial value in Eg. HI 



shock. Thus, producing the X-rays in the shocked wind of 
pulsar A in the bow shock occurring near pulsar B would re- 
quire an efficiency of ~ 0.3/C, in order to account for the 
Chand ra observa tion. 

Lvutikov ( 2004) calculated the asymptotic opening angle of 
the bow shock, a nd finds it to be 9 ~ 0.11 - 0.13 rad for the 
value of from lLvne etal] d2004). This gives Q/4n ~ (3- 
4.2) x 10~ 3 , which is in agreement with our estimate here, and 
provides an independent cross calibration for our parameter C, 
namely C ~ 2.6-3.6. 

As in ^3 the values of the hydrodynamical quantities at 
the head of the bow shock are used in order to estimate the 
emission from the shocked wind. To zeroth order, we ne- 
glect the orbital motion of the two pulsars, and their spins. 
The bow shock itself is at rest in the lab frame, in our ap- 
proximation. 7 The expressions for the hydrodynamical quan- 
tities are similar to those in |2] just that here the distance 
of the head of the bow shock from pulsar A, R = Rab - 
Rbs ~ 8.2 x 10 10 cm, is ~ 10 s times smaller. Therefore, we 
have e; nt = 4.7 erg cm J , n = 827 w5 cm J and B = lle g ' G. 
While the dynamical time in this case is much shorter, fd yn ~ 
Rbs/(c/3) = 0.6 s, the synchrotron cooling time t c oc R 2 is 
smaller by an even larger factor, so that 7 e < 7 max . Here 
7max is also constrained by radiative losses. This limit may 
be obtained by equating the cooling time t c to the accel- 
eration time, f acc = A(2Trm e cy e /eB) where A > 1, 7max,2 = 

(3e/Aa T B) l l 2 = 1.4 x W 1 A~ l l 2 e~g ,A . The limit discussed in 

13 now reads 7 maXj i = 1.2 x 10 7 ey 2 (?7/0.3), and we have 
7max = min(7 max j , 7 max , 2 ). 

Since the bow shock around pulsar B is much more com- 
pact than the bow shock on the ISM, one might expect inverse 
Compton scattering of the synchrotron photons to be more im- 
portant in this case, and therefore we check this. The Comp- 
ton y-parameter is given by Y ~ T Tl^~ l ll~ p for 2 < p < 3, 

Y - T7-7m7 r[l + ln(7m ax /7r)l for P = ^ nd Y ~ T T jP~ l J c J^Z 

for p < 2 (Panaitescu& Kumar 2000). We expect p w 2, 
for which Y(\+Y) ~ 10 _2 (r;/0.3)e e /e B . For our values of 
£(• ~ £b ~ 1, this gives Y ~ 10~ 2 , so that inverse Compton 
scattering is not very important, and is neglected in our treat- 
ment. 

We have j c = 1.1 x 10 7 e^ 1 and 7„, is still given by Eq. |2] 
The corresponding synchrotron frequencies are 

/^ m = 86g 2 4 /2 e 2 7 2 . 5 eV, (10) 
/!i/ c =17(l+F)" 2 e" ? 3/2 MeV , (11) 
/ii' max = min[30A- I (l+yr I ,20eB /2 (?7/0.3) 2 ]MeV . (12) 

The X-ray luminosity is then Lx = /xe e e r ad(^/47r)L J 4. For 

p re 2 typically e rad - 0.2, and f x - (^/^) (3_p)/2 ~ 0.02 e 3 B /4 . 
Altogether, and assuming VL/Att ~ 4 x 10~ 3 (C ~ 4), we have 

Lx ~ 10 29 e^ 4 erg s" 1 , which is a factor of ~ 20e B 3//4 smaller 
than the observed value. 

It might still be possible to increase Lx if somehow v c 
could be lowered, since this would significantly increase fx, 
and also somewhat increase e rac j. This could potentially be 
achieved if fdyn or the magnetic field experienced by the 

shocked electrons are increased (as for p pa 2, fx oc v c l ' 2 oc 

7 We ignore the slower binary period timescale, which would cause inertial 
effects, centrifugal and Coriolis, etc., as well as the possible time variability 
due to the rotation of the pulsar B magnetosphere. 
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tdynB 3 ^ 2 )- This might happen if a reasonable fraction 8 of the 
shocked wind becomes associated with closed magnetic field 
lines of pulsar B for one or more rotational periods of B 
(where Pb is ~ 4.6 times larger than the estimate we used 
for fdyn, i.e. Rb s /(c/3) « 0.6 s). In this case, this material will 
also pass through regions of higher magnetic field strength. 
This could be the case if, e.g., interchange instabilities cause 
mixing of the two fluids across the contact discontinuity. 

One might expect a modulation in the emission with the or- 
bital period due to th e change in the line of sight w.r.t. the bow 
shock ( Arons & Tavani 1993). The shocked wind is expected 
to move away from the head of the bow shock with a mildly 
relativi stic velocity, in a direction roughly parallel to the bow 
shock (Lvutikov 2004). This might cause a mild relativistic 
beaming of the radiation emitted by the shocked plas ma, re- 
sulting in a mild modulation (by a factor < 2, Arons & Tavani 
1993) of the observed emission as a function of the orbital 
phase. Another possible source of modulation with the or- 
bital period may arise if the luminosity of th e pulsar A wind 
depen ds on the angle from its rotational axis dDemorest et al.l 
2004). In this case the wind luminosity in the direction of pul- 
sar B will vary with a period P 01 -b. The duration of the Chandra 
observation, 10 4 s, is close to the orbital pe riod P m b = 2.45 hr, 
and it showed no evidence for variability (McLaughlin et al. 
2004). However, the small number of photons (77 ±9) does 
not allow to place a strong limit on a possible modulation 
with the orbital period, which might still have an amplitude 
of<50%. 

The rotation of pulsar B, assuming some misalignment of 
its magnetic pole relative to its spin axis (as expected from 
the detection of its pulses), would cause a periodic change in 
ft (with a periodicity equal to the spin period Pb = 2.77 s), 
with a n amplitude which is typically of order unity (iLvutikovl 
2004). The distance of the bow shock from pulsar A hardly 
changes, and therefore the values of the thermodynamic quan- 
tities in the shocked wind and the resulting values of fx and 
e la d should vary with a smaller amplitude. Thus, the modula- 
tion in Lx is expected to largely follow that in ft, and have a 
similar amplitude (typically of order unity). 

4. DISCUSSION 

8 In the bow shock of the solar wind around the Earth only ~ 10 -3 of 
the wind particles get captured by the Earth's magnetic field. However, the 
situation there is different in several respects from our case. For example, 
the Earth's magnetic field is nearly aligned with its rotational axis, while the 
solar wind is Newtonian (~ 400 kms~') with relatively low magnetization 
and includes protons and electrons in roughly equal numbers. Therefore this 
fraction might be larger in our case, and could possibly be sufficiently large 
for our purposes, although this is uncertain. 



Particle acceleration is expected in the binary pulsar sys- 
tem J0737-3039 both from the bow shock of the pulsar A 
wind as it interacts with the interstellar medium (ISM), and 
the bow shock of the wind of pulsar A interacting with the 
magnetosphere of pulsar B. The rotational energy loss rate, 
the systemic velocity and the orbital separation determine the 
effective angles subtended by these bow shocks, as well as the 
synchrotron peak energies in the forward and reverse shock 
systems and the radiation efficiencies at various frequencies. 
In th is model, th e likeliest explanation for the Chandra emis- 
sion (McLaughlin et al. 2004) is the pulsar A wind just behind 
the bow shock caused by the systemic motion in the ISM. In 
this case, we predict a power law spectrum which extends up 
to < 60 keV. 

The eclipse of the pulsar A radio emission near superior 
conjunction is best explained as synchrotron absorption by 
the shocked pul sar A wind in the bow shock around pulsar 
B dKaspi et al. Il2004t iDemorest et alJl2004l Lvutiko^l2004l 
lAronsetal 12001 " This explanation requires a relatively large 
number density of pairs, which in turn requires a relatively 
low wind Lorentz factor, 7,,, < 100. However, the X-ray emis- 
sion from both of the bow shocks is not very sensitive to the 
exact value of 7 W , and j w ~ 10-100 would only lower the ra- 
diative efficiency e ra d and the X-ray luminosity Lx by a factor 
of ~ 2 (for p~2) compared to j w ~ 10 5 . 

An alternative explanation for the X-ray emission, is 
simply emission from pulsar A ( McLaug hlin et al.l 120041: 
Zhang & Harding 2000). 9 In this case a large part of the X- 
ray emission is expected to be pulsed with a period Pa- In 
contrast, the emission from the bow shock around pulsar B 
might be modulated 10 at P or b or Pb, while the emission from 
the bow shock with the ISM is not expected to be modulated 
but might be angularly resolved by Chandra. 
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9 Emission from pulsar B is unlikely, since £ lot .B ~ Lx , which would re- 
quire a very high efficiency in producing X-rays. 

10 Although we find that the emission from the bow shock around pulsar 
B is likely to contribute only a few percent of the total X-ray luminosity from 
this system, it can still produce an overall modulation of up to several percent, 
which might still be detectable. 



REFERENCES 



Arons, J., Backer, D. C, Spitkovsky, A., & Kaspi, V. M. 2004, in Binary 
Radio Pulsars, ed. F. H. Rasio & I. H. Stairs (San Francisco: ASP), in press 
iastro-ph/0404159i 

Arons, J., & Tavani, M. 1993, ApJ, 403, 249 

Demorest, P., et al. 2004, submitted to ApJL 1 astro-ph/0402025 1 

Gallant, Y.A., & Arons, J., 1994, ApJ, 435, 23C 

Kaspi, V.M., Ransom, S.M., Backer, D.C., Ramachandran, R., Demorest, P., 

Arons, J. & Spitkovsky, A. 2004, submitted to ApJL 1 astro-ph/040I614 1 
Lyne, A.G., et al. 2004, Science, 303, 1 153 



Lyutikov, M. 2004, submitted to MNRAS 1 astro-ph/0403076 1 

McLaughlin, M.A., et al. 2004, ApJ, 605, L4T 

Panaitescu, a., & Kumar, P. 2000, ApJ, 543, 66 ^_^^ mp ^^. 

Ransom, S. M., et al. 2004, submitted to ApJL 1 astro-ph/0404149 1 

Spitkovsky, A., & Arons, J. 2004, ApJ, 603, 669 

Zhang, B & Harding, A.K., 2000, ApJ 532, 1150 



